On Discrete Painleve Equations 
Associated with the Lattice KdV Systems 
and the Painleve VI Equation 



F.W. Nijhoff 

, Department of Applied Mathematical Studies 

^ I The University of Leeds, Leeds LS2 9JT, UK 

■ A. Ramani 

Q 

' CPT, Ecole Polytechnique 

CNRS, UPR 14 
91128 Palaiseau, France 

B . Grammaticos 

GMPIB, Universite Paris 7 
2 Place Jussieu, Tour 24^14; 5eme etage 
75251 Paris, France 

Y. Ohta 

c/3 . Department of Applied Mathematics 



o 

(N 
00 



Faculty of Engineering, Hiroshima University 
1-4-1 Kagamiyama, Higashi-Hiroshima, 739 Japan 



1 



Abstract 



A new integrable nonautonomous nonlinear ordinary difference equation is presented 
which can be considered to be a discrete analogue of the Painleve V equation. Its derivation 
is based on the similarity reduction on the two-dimensional lattice of integrable partial dif- 
ference equations of KdV type. The new equation which is referred to as GDP (generalised 
discrete Painleve equation) contains various "discrete Painleve equations" as subcases for 
special values/limits of the parameters, some of which were already given before in the lit- 
erature. The general solution of the GDP can be expressed in terms of Painleve VI (PVI) 
transcendents. In fact, continuous PVI emerges as the equation obeyed by the solutions 
of the discrete equation in terms of the lattice parameters rather than the lattice variables 
that label the lattice sites. We show that the bilinear form of PVI is embedded naturally 
in the lattice systems leading to the GDP. Further results include the establishment of 
Backlund and Schlesinger transformations for the GDP, the corresponding isomonodromic 
deformation problem, and the self-duality of its bilinear scheme. 



Key Words: Discrete Painleve Equations; Painleve Transcendents; Similarity Reduction; 
Ordinary & Partial Difference Equations; Isomonodromic Deformation Problems; Hirota 
Bilinear Equations. 







1 Introduction 



Discrete Painleve equations ^ have now been around for already a couple of years and 
they form a new exciting class of equations, possibly leading to the definition of new 
transcendental functions. Their study may yield new insights into the analytic theory 
of ordinary difference equations, the development of which, cf. |^], has been greatly 
lagging behind compared to its continuous counterpart -the theory of ordinary differential 
equations. In view of this, the recent interest in the investigation of integrable discrete 
systems may open a new window on what might ultimately become a general theory 
of discrete and difference equations. What the study of integrable discrete systems has 
provided so far is a rich class of interesting examples of nonlinear ordinary and partial 
difference equations that are amenable to exact and rigorous approaches for their solution. 
With these examples we expect that insight can be gained into the general theory, and 
that they will help us to develop methods for the study not only of integrable equations 
but of difference equations in general. Integrable partial difference equations have been 
investigated over the last two decades, cf. [^-0], cf. |^] for a review, providing a large 
reservoir of examples for the study of further properties. To this many novel results and 
important insights have been added especially during the last few years comprising of: 
the study of initial value problems on the space-time lattice, ||9|, the study of similarity 
reductions on the lattice, the establishment of the remarkable singularity confinement 



property of integrable discrete systems, |11|, i.e. an effective integrability detector for 
discrete systems (deemed to be an initial step towards finding a purely discrete analogue 
of the Painleve property), the elaboration of proper discrete analogues of the Painleve 
transcendental equations PI-PVI, [|^, |^], the construction of special solutions and 

properties of the discrete Painleve equations, ([|| for a review). 

Let us provide a few examples. The discrete Painleve I (dPI) equation 

Xn+1+ Xn + Xn-l = — + a , = an + P , (1.1) 



arose from the theory of orthogonal polynomials (cf. ||15| for a review), and in physics in 



the matrix model approach to two-dimensional quantum gravity, [^], cf. also |13]. The 
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discrete Painleve II (dPII) equation 

I Cn^n + /I OA 

Xn+1 + Xn-1 = -I _ 2 

was found in connection with unitary matrix models of quantum gravity in and from 



the lattice similarity approach introduced in |10|. Other discrete Painleve equations have 
been found from the Backlund and Schlesinger transformations of the continuous Painleve 
equations, cf. |14|, for example the alternative discrete Painleve II (alt-dPII) equation 



Cn+l , Cn J. . 1 . /I o\ 

-\ —r=Cn-Xn-\ ha, (1.3) 



_1_ 1 _1_ 1 



that was studied more extensively in [18|. Many of these discrete Painleve equations have 



been shown to possess an associated isomonodromic deformation problem, cf. [|ig|]-|22| 



(with the notable exception of the discrete dPIV and dPV equations found in 1 12] , for which 
so far no monodromy problem has been established). Thus, we have now at our disposal 
a large body of examples of discrete equations that can be considered to be the proper 
analogues of the continuous equations. There does not yet, of course, exist a full-fledged 
theory on the discrete level, but already many of the outstanding features of discrete 
Painleve equations have been established, such as the existence of special solutions for 
special values of the parameters, the associated bilinear forms, the Backlund/ Schlesinger 
transformations, etc. (For a recent review, cf. Q). 

Central in this paper is a system of partial difference equations which can be thought 
of as a (scaling) symmetry reduction of the lattice MKdV equation, cf. |1^, 23 1. It consists 
of the following coupled set of equations: 

/ Vn,m+1 Vn+l,m \ I ^n+l,m I'n.m+l \ a \ 

P = ^ ' 

\Vn+l,m+l Vn,m ) yi'n+l.m+l ^^n,m / 

which is the lattice MKdV equation, together with 

n ■ ■ h m — ■ ■ = fi — X[—l) ^ , (l-4b) 

'^n+l,m ~l~ ^n— l,m '^n,m+l ~l~ '^n,m— 1 

which is a nonlinear nonautonomous constraint expressing the scaling symmetry. The 
similarity constraint ( 1.4b| ) was found already some years ago in ||l^ for the special values 



of the parameters A = 0, ^ = 0. Related systems were developed recently, |p3|, and a 



geometrical interpretation was given in [24|. 



It is shown that eqs. ( |l.4a ) and ( |l.4b ) are compatible on the lattice and that their 



reduction gives rise to a combined set of ordinary nonlinear nonautonomous difference 
equations of the form 

+ TT-TT = H + iy + Cn+l +Cn + 



_ ^)(^2 _ 1)^^ ^ _ ^2) ^1 (^^ ^ ^ i(-l)"(Cn - Cn+l - 2m) ^ 

^ (r + X„)(l+rX„) ' ^^-^^ 

with parameters r = q/p, m, fi and the alternating parameter v = A(— l)"""*"™, includ- 
ing the discrete-time step. We will refer to this discrete equation as GDP (generalised 
discrete Painleve equation), since as we will show in what follows this second-order non- 
linear nonautonomous difference equation incorporates many of the so far known discrete 
Painleve equations of difference typeB. In fact, special limits of (1.5) reduce to both dis- 
crete versions of PII, namely ( |1.2[) and ( |1.3D . However, by itself, eq. (1.5) can be considered 
to be a discrete analogue of the PV equation by virtue of the fact that, as we shall show 
below, eq. (1.5) reduces to PV by a continuum limit. Thus, eq. (1.5) seems to be a very 
rich equation, and this point is made even stronger by noting that solutions of (1.5) can 
actually be expressed in terms of PVI transcendents, i.e. solutions of the PVI equation 
d^u, 1 / I 1 1 \ fdwy /I , 1 1 \ dw 



+ 7 + : ^ - T + — r + 



dt"^ 2 \w w — 1 w — t J \ dt J \t t — 1 w — t J dt 

, w{w-l){w-t) ( t t-l l^ t{t-l) \ 

+ t^{t-iY r - + ^ (^;;3Tj2 - ('^ - 2^ (^^^37)2 j • (1-6) 

An exhaustive list of Schlesinger transformations for PVI has been given by Mugan and 



Sakka in |25], and in principle the GDP equation (1.5) can be considered to be the per- 



mutability condition for the corresponding Backlund transformations in the spirit of |14]. 
Also in the early paper difference relations associated with the Painleve equations 
were considered, but no explicit closed- form discrete equations associated with PVI was 
given. The GDP (1.5) is, in fact, the appropriate difference equation that is naturally 
associated with PVI and the aim of the present paper is to study this equation, as an 
integrable ordinary difference equation in its own right. 



^ There exist also discrete Painleve equations of q-difTerence type, such as the discrete PIII equation 
equation found first in jlij from the singularity confinement analysis and for which a g-difference mon- 
odromy problem was derived in Further q-difference dP's include discrete PV, [|l^ and discrete PVI 
found in S. 
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Let us point out that referring to equations like (1.5) as a discrete PV (dPV) equation 
risks to become slightly confusing. Since, on the one hand, there are various different 
discrete Painleve equations that tend to one and the same continuous Painleve equation in 
limits, on the other hand the same discrete equation may have various limits to different 
continuous and/or discrete equations, this type of nomenclature (referring to discrete 
Painleve equations in terms of the continuous PI-PVI equations to which they constitute 
the discrete analogues) is slightly deceptive. In particular, there exist various alternative 
"dPV" equations as well and so far the interconnection between these various discrete 
equations is far from clear. Issues of classification and universality will not be addressed 
in this paper, but we hope that revealing the structure behind equations such as (1.5) might 
be of help in establishing the connections between various integrable ordinary difference 
equations. 

In this paper, we will study the lattice system ( p..4| ) as well as the GDP equation (1.5) 
from various angles. The outline of the paper is as follows. In section 2 we will present 
a general formalism in terms of an infinite matrix structure from which one can derive 
the lattice equations as well as their similarity constraints. The explicit forms of the 
lattice systems are presented in section 3, together with their general solution in terms 
of PVI transcendents. In section 4 we discuss the integrability of the lattice systems 
and their consistency from the point of view of initial value problems. Furthermore, 
from the structure developed in section 2, we derive their isomonodromic deformation 
problems. Finally, in section 5 we discuss the lattice equations from a general point of 
view studying various properties, such as their Backlund and Schlesinger transformations 
as well as the general bilinear scheme underlying these equations. In fact, writing the 
lattice equations in bilinear form we are led to a system of four-dimensional discrete 
equations in which parameters and variables live on equal footing. This self-duality of 



the bilinear forms, the so-called Grand Scheme of |27] where it was observed in another 
context, seems quintessential to the integrability of the lattice system. 
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2 Derivation: Infinite-Matrix Scheme 

We start by developing a general scheme from which we obtain the discrete equations and 
their similarity reduction. The scheme is based on a structure introducing a dynamics 
in terms of infinite matrices, and it was used in the past to construct integrable partial 
difference equations together with their Lax pairs, Miura- and Backlund transformations 
and to establish the underlying algebraic structure, cf. ||6|, |^, cf. also the more recent 
account for a review. The scheme is formal, and we would like to stress that the 
objective of this derivation is exclusively to establish the basic equations and to unravel 
the underlying algebraic structure. We do not make any claim on the analytic aspects of 
the scheme: in order for it to be used to construct solutions of the equations one needs to 
establish the existence of the objects under question beyond the formal level. However, the 
scheme has proven to be highly effective in arriving at all the necessary relations between 
the equations which then in the next stage can be used for constructing analytic solutions. 

2.1 Definitions 

For brevity, let us formulate the structure in a rather abstract way. We basically need 
three ingredients to formulate the entire scheme: 

• An infinite ( Z x Z) matrix C which we can take of the form 



depends on additional variables that are to be determined later. The integrations 
over contour T and measure dX need not be specified at this point but we we will 
loosely assume that they can be chosen such that the objects to be introduced below 
are well-defined. From the definition it is clear that C is symmetric: C = *C , (the 
left superscript t denoting the adjoint of the Z x Z matrix). 

• Matrices A and *A that define the operations of index-raising when multiplied at 
the left respectively the right, as well as matrices I and 'I that count the index 





in which C£ and 



are infinite vectors with components {ci)j = (*C£)j = P , and p£ 
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label, namely by (I • Ci)j = j(ci)j, (*c^ • = j{*Ci)j. The interrelation between 
the variables A, *A,I and is given by 

A^' • I = (I + j) ■ , (2.2a) 
*I • W = *A • ( *I + j) . (2.2b) 

• An infinite matrix fl obeying the equations: 

ftA^ - (- *A)Jn = Oj , (2.3a) 
*I-n + n- I + n = 0. (2.3b) 

in which 

fe— 1 
j=0 

and where O is a projection matrix on the central element, i.e. (O • Ck)i,j = '^i,oCoj, 
etc. , obeying also 

*I • O = O • I = . (2.5) 

There are linear equations for the object C, the form of which depends on the choice 
of the factor p£. We can choose pi to depend on continous time-variables ti, ^3) • • • j in 
terms of which we have linear evolutionary flows. However, since we deal here with the 
discrete case we choose pn of the formi 

Thus we can impose for C the following type of discrete evolution 

C • (p - *A) = (p + A) • C , (2.7a) 
C-((7-*A) = (g + A).C, (2.7b) 

in which the shift C ^ C is the shift in the discrete variable n associated with a 
(complex) parameter p and the shift C ^ C is the shift in the discrete variable m 



■^Throughout this paper we wiU use independent discrete variables n^, or in particular n and m. Strictly 
speaking, these variables need not be integers themselves, we only require that they shift by integers, i.e. 
riv (z 9^ + Z, where the 9i, are arbitrary constants. 
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associated with the (complex) parameter q. In principle we can select an arbitrary set 
of lattice parameters p^, each associated with its own lattice shift, and in view of the 
linearity of the corresponding equations for C we can impose all the discrete evolutions 
for all chosen values pi, simultaneously. The true objects of interest, however, will obey 
nonlinear equations, and these objects are the following: 

• The infinite matrix 

U = C-(l + rj-C)"^ , (2.8) 

• An infinite determinant 

r = detz (1 + • C) . (2.9) 
To make sense of the latter we can use the expansion 

det(l + ^) = l + ^Aii + Y^ 

i i<j 

and imposing that the integrations in ( |2.1| ) are such that the terms in the expansion which 
are all of the form 

truncates. The relevant quantities in terms of which one can derive closed- form nonlinear 
equations are actually the individual entries of the infinite matrix U, which -by the way- 
is symmetric, *U = U as a consequence of the symmetry of C. 

2.2 Basic Relations 

To derive the basic equations from this scheme, we mention the following relations that 
can be derived from the definitions (for details, cf. |I]). For the r- function we have 

J = det z (l + (P - *A)'' • O • U) = 1 + (U • (p - *A)~')^^g , (2.10) 

whereas the infinite matrix U obeys the matrix Riccati equation 

fj. (p- *A) = (p + A) -U-U- O -U . (2.11) 



A ■■ A ■ ■ 
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Eq. ( 2.1l| ) forms the starting point for the construction of integrable lattice equations. By 



combining different shifts associated with different parameters p one can actually derive all 
relevant discrete equations within the KdV family from this single relation. Eq. ( |2.11| ) has 
its continuous counterpart as well which can be derived by a similar set of manipulations 
starting from 



Once we have selected a particular set of continuous or discrete variables we have yet 
another type of relations that we can consistently impose on the quantity C, involving the 
operation generated by the matrices I and *I. 

The dependence of the plane wave-factor pi on both the discrete variables riiy as well 
as on the lattice parameters piy and the spectral variable i is what determines the various 
equations for the relevant quantities. Clearly, from the definition (|2.6|) these are given by 



nPi = ] Pi , (2.12) 



for the dependence of pe on the n,^, (the Ti, denoting the shift with respect to the variable 
TT-jy, i.e. Tupi{ny) = p^{ny + 1)). Furthermore, the dependence on the lattice variables is 
given by 

^ Pi = nu ( — 3_ ^ (2.13) 



whereas the dependence on the spectral variable is given by 



4^^ = E-^P^(^-^)/'^- (2-14) 



Note that from ( |2.12| ) and (2.13) we immediately have that the operations of shift in 



n 



V 



and differentiation w.r.t. p^ commute: 

Keeping in mind the representation ( p.l| ) for the infinite matrix C, we are led to the 
following choice of linear constraint: 



= C + IC + C*I 

ar 

+ Y^p,nJc- ^-C) , (2.15) 



reflecting the scaling invariance. The term on the left-hand side is the boundary term from 
the integrations over the curve T in the complex ^-plane. We don't want to assume on this 



formal level too much concerning the analytic details of the functions involved within the 
integrations, and thus we will make the assumption that the boundary terms are primarily 
dominated by the behaviour of the integrants in the limits as i tends either to zero or to 
infinity. Thus, formally, we assume that the boundary terms can be interpreted in the 
following way 

[/(^)]ar = f A hm -/xlim)/(^) , 

for any "regular" function /. It will turn out that with this assumption these terms will 
lead precisely to the correct behaviour accounting for the appearence of the free parameters 
in the corresponding nonlinear equations as a consequence of the asymptotic behaviour of 
the factors pi as £ ^ and as £ — > oo, namely 



l + 2ej:^{n^/p^) + 0{f) as 1^0 



In addition to (2.15) we have the dependence on the lattice parameters, in terms of 
which C obeys the following linear differential equations 



as a consequence of which the discrete equation will be compatible with continuous equa- 
tions in terms of the lattice parameters. To derive the corresponding equation for U we 
need to multiply (2.15) from the left by (1 + • C)"^ = 1 - • U and use eq. (U) to 
derive the following nonlinear counterpart of (2.15) 



= U + I- U + U-*I 

ar 

+ X! "^^P^ - *A Pu+ ^. ^ ^ Py - *A ^ Py + K ^ , 

(2.17) 



as well as the nonlinear counterpart of (2.16) 



^U = n, f^--U-U- ^ + U- L^.O-^— --u) . (2.18) 

opv \Pu + A Pu - ^ Pv- rs. P;. + A ) 

In the boundary terms of eq. (2.17) we have introduced the infinite vector 

= (1 - U • f)) • . (2.19) 



It will turn out that the vector will be the starting point for the construction of the 
monodromy problems for the resulting lattice systems. In fact, components of this vector 
will form the eigenfunctions of the corresponding Lax pairs, and the variable i will play 
the role of the spectral parameter. 

Now that we have the framework in place from which the similarity reduction on 
the lattice can be derived, we can exploit the relations derived earlier to get closed-form 
discrete equations, as well as their similarity reduction. The KdV family is characterized 
by the symmetry 

C = *C U = *U , 

leading to an additional set of algebraic relations of the form 

U • (- *A)^' = A-'' • U - U • Oj • U , 3 even . (2.20) 

This condition on the infinite matrix system will be specific to the KdV class of lattice 
equations. In the next sections we will employ the various relations for the infinite matrix 
U to derive closed-form equations for some of its entries leading to on the one hand 
integrable lattice equations as on the other hand their similarity constraints in closed 
form. 



2.3 Lattice Equations 

Let us now introduce several objects in terms of which one can derive from the basic 
system ( p. 11 ) closed- form equations, namely 



u.U„,„ , v.-|^-U^^)^^, (2.21a) 



as well as 



v^ = l- -U] , Sa = a- • U • *a) . (2.21b) 

Va + A /o,o Va + A /o,o 

In ( |2.21| ) a and (3 are arbitrary parameters which we can choose at our convenience. The 
KdV case is distinguished by the fact that the variable Sa,/3 is symmetric in the parameters 
a and /3, i.e. Sq,^^ = sp^a- It is obvious from the definitions that Va can be obtained from 



10 



Sa,i3 in the limit (3 ^ oo , and that u is recovered in the simultaneous limit a, f3 ^ oo . 
However, we shall be interested as well in cases where either q = or /? = or both. 



From the basic equations ( |2.11| ) we can now derive the following set of relations for 
the objects defined in ( ^.21 ): 



as well as 



Sa = {P + u)Va - ip- a)Va 

•S/3 = (p + (3)vp -{p- u)vp , 

1 - (p + P)Sa,f3 + {p- a)Sa,l3 = VaVp 



(2.22a) 
(2.22b) 

(2.23) 



in which the tilde denotes the shift Ty with respect to a discrete variable riy = n associ- 
ated with the lattice parameter py = p. Introducing a second discrete variable = m 
associated with the lattice parameter pi, = q, we have two different transformations on 
the lattice, namely 

Pk^Pk , Pk'-^Pk, 

which can be combined in order to obtained closed-form equations. Since the dependence 
on the two variables n and m enter in C via the pi of ( |2.6| ), it is clear that C depends 
linearly on these variables implying the permutability of the lattice shifts. This allows us 
by combining two lattice shifts to derive closed- from partial difference equations. In fact, 
all objects ( |2.21a ) and ( 2.21b| ) are then interpreted as functions of the lattice sites (n,m), 
and the shifts indicated by ~ and ^ represent the translations 

(n, m) 1-^ (n + 1, m) , (n, m) i-^ (n, m + 1) . 

Thus, for the objects (|2.21a|) and (|2.21b| ) we derive from the various relations (|2.22| ) and 
( |2.23| ) nonlinear partial difference equations leading to the various lattice equations within 
the KdV family. 

The most general lattice equation that can be derived directly from ( |2.23| ) is 

l-{p + f3)sa^i3 + {p- a)sa,f3 ^ 1 - (g + a)sa^i3 + {q- P)Sa,l3 ^4) 

l-{q + P)Sa,l3 + {q- a)Sa,f3 1 - (P + a)^a,l3 + {P - /?)««,/? ' 

cf . 1^] , which for fixed values of a and (3 is an integrable partial difference equation for Sa,i3 
with as independence variables the lattice sites {n,m). The p and q can be interpreted as 
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the corresponding lattice parameters. It must be pointed out that the independent discrete 
variables n and m, are attached to the lattice parameters p respectively q, so that fixing the 
latter means choosing certain lattice directions of an in principle multidimensional grid. 
As is clear from ( |2.6D , we are absolutely free to introduce additional lattice directions each 
carrying its own lattice parameter. In terms of each pair of these lattice directions we have 
an equation of the form ( p. 24 ) for one and the same object Sa,i3, and all these equations can 



be shown to be compatible. From that perspective the equation ( 2.24 ) actually represents 
a parameter-family of compatible equations, which is the precise analogue of the hierarchy 
of nonlinear evolution equations in the case of the continuous KdV systems. 



Eq. ( 2.24| ) has several specialisations according to the choice of the latter parameters. 



The equations satisfied by the variables u and Va follow from the limits a ^ oo, /3 — > oo , 
respectively a arbitrary, /? — > oo . Between these variables there are relations of Miura 



type which can be derived from ( 2.22 ), eliminating the intermediate variable Sq, by using 



both translation on the lattice, namely 



p — q + u — u = {p — a):^ — {q — a)^ 

Va Va 



{p + P)^-{q + P)^ , (2.25a) 

V/3 Vp 



p + q + u — u = ip — a)-:^ + {q + a)^ 

Va Va 

= {p + P)^ + iq-P)'^ ■ (2.25b) 

V/B Vp 

Within the lattice KdV family we will concentrate on the three main equations, namely 
the ones for for u, vq and z, where z is related to so,o by 

^ = so,o----, (2.26) 
p q 

because it is for these variables that closed-form similarity constraints can be derived. The 
lattice equations satisfied by these objects can be easily deducted from the relations given 



above or directly from ( 2.24 ) by the appropriate limits. The explicit forms will be given 
in section 3.1. 



From the point of view of considering (2.24) and its atavars as parameter-families of 



compatible lattice equations, it makes sense also to consider the continuous equations 
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expressing the dependence of the variables u, v and s on the lattice parameters. The 
resulting differential equations are by construction compatible with the lattice equations 
and integrable by themselves. They follow from the basic equations 
du 

— = n (1 - VpV-p) , (2.27a) 

— — = -n h Sa-pVp , (2.27b) 

op \ p — a ) 

^ = n f ^^^^^^ - £-/l^£p^ - s^^.psp^ , (2.27c) 
op \ p — a P + P / 

leading in particular to the differential-difference equations for the variables with a,j3 = 

0,00, namely 

?^ = nfl-;^^l, (2.28a) 
dp \ 2p + u-uj ^ ' 

— log uo = ^ , , (2.28b) 

op PVQ^ Vq 

|£ = _2n (?-.)(.-.) ^ 

Obviously, similar equations hold in terms of the parameter q. 

Finally, we mention the bilinear equations that are obeyed by the r- function, that 
was formally defined by the representation (12.91) from which one can deduce the following 



identifications of the variables Vp and Vq in terms of the r-function 

Vp=^^ , vq=^^. (2.29) 

(The undertilde ^ and underhat denote the lattice translations in the reverse directions to 
the shift ~respectively^). In fact, we can associate with every value of the parameter a, in- 
stead of p or g, a variable Va by similar expressions as ( |2.29| ) in terms of lattice translations 
associated with that value. In particular, with a = we obtain the identification 

t^o = — , (2.30) 

in which 

r± = det z (1 ± f7 • C) . (2.31) 

The r_ can be viewed as the r-function obtained from r = r+, i.e. ( |2.9D , by applying a 
lattice shift associated with parameter value a = 0. For both t± we have the following 
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bilinear lattice equations 



{p + q)t'Z = 2g'TT + {p - q)TT , 
ip + l)^Z — 2prf — {p — q)z^ ■, 



(2.32a) 
(2.32b) 



which, together with the identifications ( 2.29| ), are easily obtained from the equation for 
Vp, namely 

(2.33) 



{p + q)J^ + {p- q)^ = '^p , 



respectively from a similar equation for Vq. The bilinear equations (2.32a) and (2.32b) are 
not independent but there is a consistency that follows from the following diagram 







1 ( 


i K 

I §: 



Figure 1: Consistency of the bilinear equations. 

In fact, taking initial data at the vertices of the first diagram in Fig. 1, we can calculate 
the remaining vertices in the rectangle vertically by using (2.32a) and horizontally by using 
(2.32b). Then the value of r at the vertex indicated by (CS) at the right-bottom corner 
can be calculated in two independent ways, but it is easily checked by explicit calculation 
that this does not lead to an inconsistency: we obtain the same value irrespective of the 
way in which we calculate that value. 

Mixed equations between t+ and can be derived from formulae like: 

^^ = (P + .)I , P-^ = {p-qf^, (2.34) 



which follow from (2.25a) and (2.25b) for different values of a and f3. Eqs. ( 2.34| ) lead to 
the following bilinear relations 



{P + q)r±T^ = PT±T^ + qT^T± . 
{P - Q)t±T^ = PT±T^ - qT^T± . 



(2.35a) 
(2.35b) 
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For T we have also differential relations in terms of the lattice parameters. In fact, 
from the definition ( |2.9| ) one can easily deduce that 

d 

— logT = nSp-p, (2.36) 
and from ( 2.23| ) for a = —p, (3 = p, together with ( |2.29| ) we then have the bilinear equations 



of the form 



X ^ 2p^dT 2p dr 

TT = TT-\ T- tt'^^ ' (2.37a) 

n op n + 1 op 

TT = TT-\ r— I ) (2.37b) 



n \ dp dp 

which will provide us with bilinear equations of differential-difference type. 
2.4 Similarity Constraints 

To obtain closed- form constraints for single objects we need the cases that the parameters 
a and [3 go either to zero or infinity (otherwise we will get constraints containing derivatives 
with respect to the variables a and (3). So, we will focus primarily on similarity constraints 
for the objects: u, vq and sq.o, but we need some other objects along the way as well. The 
similarity constraints for u, vq and z are obtained from 

2A(-l)"+'"(np + mq-u)=u + np{vpV_p - 1) + mq{vqV^g - 1) , (2.38a) 
{fj. - A(-l)"+™) vo = n [(1 + pso-p)vp -vo]+m [(1 + qso-q)vg - vq] , 

(2.38b) 

Tl Tfi 

-(2/i + l)z = -(1 + pso -p)(l - pspfl) H {1 + g'So,--g)(l - gsq,o) , (2.38c) 

p q 



with ( 2.26| ), and in which we need to use the relations 



1 +psf)-p = VQV-p , (2.39a) 

1 - PSpfi = VpVo , (2.39b) 

-{p- u)v-p = s-p , (2.39c) 

and various other expressions that follow directly from ( |2.22| ), ( |2.23 ), such as 



v„V-„ = 1 , 2p + u-u = = 2p^ , (2.40) 

vo Vp 
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(and similar equations with • replaced by ^ and p replaced by q). In the next section we 
will collect the explicit forms of the resulting equations. 



situation where there are only two discrete variables n and m associated with lattice 
paramters p and q respectively. However, we are absolutely free to extend the equations 
to include other lattice variables associated with additional directions (each carrying its 
own lattice parameter) . For each of such lattice direction we have to add on a term of the 
same form as the ones present on the right-hand side of eqs. (2.38a), (2.38b) and (2.38c). 
In the physics of two-dimensional quantum gravity this has its analogue in the famous 
string equation, where it has occurred in the literature on random matrix models, cf. e.g. 
| p!3| , p!7| . Effectively, the procedure of increasing the number of variables in the similarity 
constraints amounts to building what is the discrete analogue of the Painleve hierarchy, 
which in the continuum leads to increasingly higher-order equations after the elimination 
of redundant variables. This point of view on discrete Painleve hierarchies is interesting, 
but to develop it furher is beyond the scope of the present paper. 

3 The Lattice Painleve System 

On the basis of the structure exposed in the previous section, we will present in this section 
the closed- form expressions for the similarity reductions of the members of the lattice KdV 
family, which comprises of the three classes of equations: those for the variables u, v and 
z which will be identified with respectively the lattice potential KdV equation, the lattice 
potential MKdV equation and the lattice Schwarzian KdV (SKdV) equation. On the 
continuous level these equations are given respectively by: 



Remark 2.1: In the similarity constraints ( p. 38 ) we have restricted ourselves to the 



+ 3ul 



'X 1 



(3.1a) 



- 3 



Vx'^xx 



(3.1b) 



— Vxxx 



V 
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Eq. (3.1b) is the potential MKdV equation, i.e. the variable v = {logv)x obeys the MKdV 
equation. Eq. (3.1c), (where S{z) denotes the Schwarzian derivative of z), is the SKdV 
equation which is invariant under Mobius-transformations: 

z ^ , ad-bc^O. (3.2) 

cz + d ^ ^ 

As is well-known, these equations are interrelated via Miura transformations and effectively 
form the Miura chain: 

SKdV — > MKdV — > KdV 

In fact, in the first step the connections is given by the Hopf-Cole transformation 

9 — 1 ^XX — —9 1 '^t / \ 

Zx = V =^ V = <^ Vx — V = , (3.3) 

2 Zx 2zx 

whereas in the second step we have the Miura transformation 

Ux = ls{z) = '^-2'^=Vx-v' . (3.4) 

The much-heralded conjecture of [28| postulating the connection between similarity 



reductions of integrable nonlinear evolution equations and ODE's having the Painleve 



property has been the motivation in [10| for pursuing a similar connection on the discrete 
level. Rather than trying to find directly a similarity variable (which on the discrete level is 
rather awkward) it was found to be more convenient to formulate the reduction in terms 
of constraints on the solution of the lattice equations. To illustrate the idea by means 



of the continuous equations (3J), it is clear that the similarity reduction under scaling 



symmetry of the SKdV equation, which, as was noted in [29|, turns out to be of the form 



z{x,t) = tf^Z{0 , ^ = xr^/2, (3.5) 

where ft = |(1 + 2fi) is an arbitrary constant (the parameter fi we introduce for later 
purpose), can be cast in the form of a linear nonautonomous constraint, namely 

3flz = xzx + 3tzt . (3-6) 
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Imposing ( |3.6D on the solutions of (3.1c) and eliminating the t-derivatives we obtain the 
following third-order ODE 

Eq. ( |3.7| ), in which t is now merely playing the role of a parameter of the equation, was 
referred to in |29] as "Schwarzian PII" (SPII), because of its relation to the Painleve II 



(PII) equation. In fact, via the Hopf-Cole transformation ( |3.3D we obtain from ( |3.7D the 
Painleve II (PII) equation 

v., = 27)3 - ^t; + £ , (3.8) 
which by itself can be obtained also from a linear constraint, namely 

fiv = xvx + Stvt =^ V + xvx + Stvt = , (3-9) 



which represents the well-known reduction of the MKdV to PII, cf. [Q. To complete the 
picture, we mention that the similarity constraint for the potential KdV is given by 

u + xux + 3tut = f{t) , (3.10) 

with the freedom to have an arbitrary function f{t) on the right-hand side. Thus, we 
obtain for U = Ux + fi Painleve 34 (P34), i.e. the equation no. 34 in Painleve's list, 
1 31, 52l, namely 



UUxx - -ul + 2C/3 - — + ii- = . (3.11) 
2 ^ 3t 8*2 ^ ' 

From the Miura transformation (^^) between the MKdV and KdV equation one easily 

obtains the Miura and inverse Miura transformations relating PII to P34, [^] , namely 

In the next subsection we will give the precise discrete analogues of these formulae and 
demonstrate that the entire Miura scheme, as well as its similarity reduction, survives on 
the discrete level. The main difference between the continuous situation and the discrete 
one lies in the fact that all similarity constraints become nonlinear on the discrete level, 
and therefore are not easily solved by extracting an explicit similarity variable. The actual 
forms of the constraints are highly nontrivial and follow from the infinite matrix structure 
of the previous section. Nonetheless, as we will show explicitely, all the interrelations 
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between the various discrete equations and similarity constraints will exist on the discrete 
level and can, therefore, be checked independently from the equations themselves without 
reference to the structure developed in section 2. 



3.1 Explicit Forms 

The first example of a lattice equation we consider is the cross-ratio equation | 

i^n,m •2'n,m+l) (^n,+l,m •2n+l,m+l) p'^ 

which is obtained from ( p. 241 ) by taking a = f3 = 0. This equation expresses the condition 
that under the map (n, m) Zn^m the canonical cross-ratio of four adjacent points on the 
lattice is constant, which was interpreted recently as a notion of discrete holomorphicity, 



(3.13) 



1 34 1 in connection with its invariance under ( |3.2| ). Since it is well-known that the cross- 
ratio forms the discrete analogue of the Schwarzian derivative, eq. ( p.l3 ) is the natural 
discrete version of the Schwarzian KdV equation (3.1c). 



In the discrete case the similarity constraint for the lattice Schwarzian KdV ( p. 13 ) 
takes on the form 



Zn+l,m Zji—l^m ^n,m+l ^n,m— 1 

(3.14) 



which was first given in |29[|. As was pointed out in [^], this nonautonomous nonlinear 
constraint coincides precisely with the condition that the map (n, m) i-^ Zn^m can be 
identified with the discrete analogue of the power map C C"" sls a discrete conformal 
map in the sense of [^]. This is clear by comparison with the continuum situation, where 
it resides in the fact that the time-dependent prefactor in the solution of the similarity 
constraint, i.e. (|3.5| ), is given exactly by the power map t tP" in the complex t-plane. 

Eq. ( ^.131) can be trivially resolved by identifying the differences in the variables z in 
terms of products of variables v by means of the following relations 

p{,Zn,m ^n+l,m) — Vn+l,inVn,m i (3.15a) 
Q{Zn,m •2n,m+l) — 'Vn,m+l'^n,m ; (3.15b) 
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which in turn by ehminating the z variables in an additive way leads to the following 
nonlinear equation for the new variable v = Vn,m, (which coincides with the variable vq in 
section 2), 

PVn,m+lVn,m + QVn+l,m+lVn,rn+l = QVn+l,mVn,m + PVn+l,m+lVn+l,m , (3.16) 

which is the lattice version of the potential MKdV equation (3.1b). The similarity con- 
straint (|3.14| ) can then be used to obtain a similarity constraint in terms of v. In fact, 



taking the difference of ( |3.14| ) in the n- variable and using ( |3.15| ) we obtain 

, 1 / , 1 \ ^n+2,m '^n— 1,171 

/i + - = (n + 1) ■ n 



P I '^n+l,m+lVn+l,m,—l '^n,m+l'^n,m—l 



+m- 



Q \ f^n,m(^n+l,m+l ~l~ '^n+l,™— l) "^^n+ljm (f^n,m+l ~l~ ^n,m— l) / 

which by using next ( p. 161) to transform the last two terms can be cast into the form 
H= -{n+ 1)- — + -n 



2 Vn-\-2,m ~\~ ^71, m 2 fjj-j-i^m -|- f^— l,r7i 

Vn+l,m+lVn,m+l ''^n+l,m— l^n,m— 1 

7 V V ' 

{^71+1,7X1+1 ~\~ '^71+1,7X1— l)\yn,7n+l ~\~ Vn,7n—l) 

where the right-hand side is actually the sum of two terms plus their n-shifted counter- 
parts. Thus, after performing one discrete "integration" we obtain the following similarity 
constraint in terms of fn,jn 

^ Vn+l,m - Vn-l,ra ^ ^ Vn,m+1 - Vn,m-1 ^ ^ _ ^(_i)n+m ^ (g^;^^) 
Vn+l,m ~\~ Vfi—l.m Vn.xn+l ~l~ ^n,m— 1 

where A appears as an integration constant. (Obviously, the whole argument goes equally 
through with m-shifts as well as with the n-shifts). The alternating term arising from 
a discrete integration is typical for the discrete situation and adds effectively an extra 
parameter to the equation. The resulting similarity constraint for Vn^m is exactly the one 
that one obtains from the formalism of section 2, namely eq. (2.38b). 

It is on this level of the equation for Vn,7n that one can actually eliminate the depen- 
dence on one or the other of the discrete variables n or m in the discrete case. In fact, 
concentrating on the n-dependence, we define as dependent variables 



,yn = Xn{m) = '^"+^'"' ^ (3.18) 



Vn+l,m+l Vn,7n+1 
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which through the lattice MKdV equation ( |3.16| ) are related via 



^ yn = — , (3.19) 



1 — ryn 1 + rx. 

with r = q/p, as well as the combinations 

an = an{m) = '■ — , 6„ = = — '■ — '■ . (3.20) 

Vn+l,m ~r I'n— l,m ^n,m+l ~r Vn,m—1 

In order to eliminate the m-shifts from the equations we need to express 6„ in terms of x„, 
or (equivalently) ?/„. To achieve that we add to both sides of eq. ( 3.1(^ ) the corresponding 



terms of its m-backshifted counterpart which leads to the following relation 

{Xn + r)bn+l +Xn = {Vn " r)hn + Vn ■ (3.21) 

Noting further that a„ can be easily expressed in terms of x„ (or y„) through the relation 
Vn/xn-i = Vn+i/vn-1 , leading to 

Un ~ Xn-1 „„N 
an = — , [6.22) 

Un ~r Xn—1 

we can now use the similarity constraint ( p. 17 ) to eliminate 6„ from ( 3.21| ) to obtain 



finally a closed-form ordinary difference equation in terms of either the variable Xn or y„ 
(or both), which reads 

, Xn+l - Xn + r{l - XnXn+l) 



(n + l)(r + x„)(l + rxj 
— n(l — r'^)x 



Xn+l +Xn + r(l + 

XnXn+1 ) 

Xn Xn—l ~\~ ?"(1 XnXn—l) 



Xn ~\~ Xn—1 ~\~ r(l -\- XnXn—l^ 

= ^r(l + 2rxn + xl) + A(-l)"+"*(r + 2x„ + rx^) - mr(l - xl) , (3.23) 

in which r = q/p, m, fi and A are parameters of the equation. In subsection 3.2 we show 
that (3.23) can, in fact, be regarded to be a discrete analogue of both the PHI as well as 
PV, and that in special limits of the parameters it contains both the discrete Painleve II 
equation (dPII), eq. (|L^), as well as its alternate version (alt-dPII), eq. (|L^). We conclude 
that eq. (3.23) seems to unify the various branches of discrete Painleve equations. To 
obtain the equation in the form (1.5) we proceed slightly differently, namely by redefining 
the variables according to the odd/even lattice sites and setting 

-^2n = X2n , X2n+1 = ■ (3.24) 
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Using then the relations ( ^.21 ) and ( |3.22|) , as well as (|3.19| ) for the even and odd sites, eq. 
(1.5) follow after a straightforward calculation. 

Arriving now at the equations for the variable u, we have the following partial difference 
equation 



{p- q + Un,m+1 - Un+l,m) {P + Q + Un,m - Wn+l,m+l) = P 



2 „2 



(3.25) 



which is the lattice (potential) KdV equation, which can easily derived from the relations 
in section 2. Eq. ( |3.25D resolves trivially the lattice MKdV equation, i.e. the lattice 



equation for Vn,m in the form ( 1.4a ), via the Miura relations 



p - q + Un,m+l - Un+l,m = P 



P + q + Un,m - Un+l,m+l = P- 



+ q 



- q 

Vn+l,m+l 



The similarity reduction of the KdV is obtained from the constraint 



(3.26a) 
(3.26b) 



(A(-l)"+™ + \){u„ 



np — mq) + 



np 



2p + Un-l,r. 



+ 



mq 



Iq -\- Un,m—1 '^n,m+l 



0, 
(3.27) 



noting that the free parameter. A, now enters purely in the alternating term. Thus, in 
a sense the similarity constraint for fn,m interpolates between two extreme cases: the 
equation for Zn,m in which only the parameter /i enters, and the equation for Un,m in 
which only the parameter A is present. We have already shown that the constraints ( ^.14 ) 
and (3.17) are related through the lattice Hopf-Cole transformations ( p.l5| ). It can also 



be shown that the constraints (3.17) and (3.27) are related, namely via the lattice Miura 
transformations ( 3.26| ). In fact, using relations of the type ( ^.40D we can rewrite ( |3.27| ) in 
the form 



1 



(A(-l)"+™ + ){un,m -np- mq) + 



npVn,m 



+ 



mqvn,m 



, (3.28) 



and then by subtracting from the n-shifted counterpart of ( |3.27] ) its m-shifted counterpart 
and using (3.26a) together with the lattice MKdV equation ( p.l6| ), the relation 

PVn,m + qVn+l,m+l _ ^ _^ PVn+l,m. " qVn,m+l 



Vn+l,m+l + ^^n+l,m-l 
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and a similar formula with p and q and their respective shifts interchanged, we can derive 
the following 

) = qVn,m+l {Cn ,m+l Cn,m) ; 

where Cn,m stands for the constraint, i.e. Cn.m = nan + rnbn ~ A* + A(— l)""*"'" . Alterna- 
tively, we can subtract from ( |3.28|) its n- and m- shifted counterpart and use (3.26b) to 
derive 

PVn,m {Cn+l,m ~ Cn,m) = QVn+l,m+l {Cn+1 ,m+l Cn+l^'m) ? 

and obviously there is a similar relation with p and q and the n- and m-shifts interchanged. 
Clearly, if Cn,m is identically zero all these relations are trivially satisfied demonstrating 
the consistency of the constraint ( 3.27| ) with the constraint ( 3.17 ). 



Remark 3.1: As is evident from the comparison with the continuous situation (see 
previous subsection) the system consisting of ( p. 25 ) and ( p. 27 ) can be viewed as forming 



a discretisation of the discrete P34, ( 3.11| ). The elimination of the m-shifted objects in 



favour of n-shifted ones yields an algebraic equation, and we will not give it here. 

Finally, we mention the bilinear forms for the lattice KdV family. They are given 
in terms of coupled system for two r-functions, r and a, which are identified with r+ 



respectively t_ given in (|2.3lD , i.e. T_j_ — Tn^m and t_ — o'n,m- For each these we have the 



following bilinear lattice equations 

{p + q)Tn,m+lTn+l,m-l = '2.qTn,mTn+l,m + {P - q)rn,m-lTn+l,m+l , (3.29a) 
{P + q)Tn+l,mTn-l,m+l = '^P'Tn,m.Tn,m+l - {P - q)Tn-l,mTn+l,m+l , (3.29b) 

and exactly the same equations for an,m- The coupling between r and a takes place via 
the additional relations 

(^n+l,mTn—l,m ~^ '^n—l,m'^n+l,m — '^^n,niTn,m i (3.30a) 
'^n,m+l'Tn,m—l ~\~ '^n,m—l'^n,m+l — '^^n,mTn,m i (3.30b) 

which follow from ( |3.29| ) by taking instead of the n- or m-shift a shift in a third di- 
rection associated with a lattice parameter equal to zero, which amounts to a Backlund 
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transformation Tn-m ^ crn,m- More generally, we obtain from (|2.35| ) the coupled system 

P(^n,mTn+l,m+l + q(^n+l,m+lTn,m = {P + Q)o'n+l,mTn,m+l , (3.31a) 
P(^n+l,mTn,m+l — QO'n,m+lTn+l,m = {P — Q)o'n,mTn+l,m+l , (3.31b) 

as well as two similar equations with a and r interchanged. Finally, the similarity con- 
straints, which for the r-functions are given by the coupled system of bilinear equations: 



{n + m + IX- A(-l)"'+™) <Tn,mTn,m = nan+l,m.Tn-l,m + man,m+lTn,m-l , 

(3.32a) 

{n + m- fl + A(-l)"'+™) an,mTn,m = nan-l,mTn+l,m + mcT„,m-ir„,m+l , 

(3.32b) 



which in conjunction with the lattice bilinear equations (2.32) and the relations ( p. 311 ) 



yield the bilinearisation of the lattice Painleve system, and thereby of the GDP (1.5). 
The bilinear form ( 3.32| ) of the similarity constraint of [^] were first given in |^5|. We 



return to bilinear forms in section 5 within the context of Backlund- and Schlesinger 
transformations for the general lattice Painleve system. We will there extend the bilinear 
formulation treating both independent variables n, m and free parameters A, fi on the 
same footing, and demonstrate the important phenomenon of self-duality for the system 
in this extended bilinear form. 



Remark 3.2: It is remarkable that the similarity constraint for the variable z, i.e. ( p. 14]) 
can be reformulated in various ways. In fact, since we have shown that together with the 
lattice equation ( p. 13 ) one can actually derive the similarity constraint for v, leading 
to (3.17), via the Miura relations (|3.15| ) which resolve the equation (3.13), we can then 
retransform the constraint for v, again by using the Miura relations, into a seemingly 
different similarity constraint for z, namely 



n+m 



n- 



•2n+l,m, ~l" ^n—1, 



2^ri 



Zn+l,ri 



Zn—l,m 



+ m- 



1+1 ~l~ ^n,m—l '^^n,'. 
^n,m+l ^n,m—l 



(3.33) 
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and pushing this even further, we can then combine ( |3.33| ) with the original constraint 



(3.14) to obtain the following third form of the similarity constraint 



;l + 2A(-l)"^ 



2n — ■ h 2m — ■ . (3.34) 



^n— 1,' 



Obviously, the three constraints ( 3.14 ), (|3.33|) and ( 3.34 ) are not independent since we 
have indicated how to derive one from the other. However, it remains an amusing fact 
that these three forms, each involving the Painleve parameters jj, and A in a different way, 
exist and are equally valid. 



3.2 Continuum Limits 

The general discrete Painleve equation (3.23) is a quite rich equation in that for special 
values of the parameters it contains previously known discrete Painleve equations, namely 
the discrete analogues of the PII equation!. In fact, taking the limit r — > oo, m/r ^, it 
is easy to see that (3.23) reduces to 

n + 1 n 1 1 / I \ 

H 7 = o +/^+ oM ' (3-35) 



^n^n+l ~1~ 1 •^U'^n—l ~1~ 1 2 2 \ Xfi 

which is the the alternate dPII equation (1.3). 



The limit to the dPII equation ( [1.2| ) is a bit more subtle, since we need to work in an 
oblique direction. In fact, to obtain dPII from (3.23) we write r = 1 + 6 and take the limit 
6^0, whilst taking n = n' — m, where n, m — > co such that n' is fixed (in the limit) and 
5m — > r] finite. In that limit we have also (1 — x„)/(l + x„) — > a„+i and (3.23) reduces to 
the following equation: 

in' + 1)(1 - xl) + {-Xr')il + xnf = 27?x„ (i^^ + 1^^) , (3.36) 

where we consider x„ to be a function of n' rather than n. In terms of a„, (omitting from 
now on the primes on the n' variable), this equation reads 



1 / , N -/X + A(-l)" + na„ 
77??(an+i + an-i) = 1 5 , (3.37) 

2 1 — 



^For this reason we have often in private referred to (3.23) as the master d-Pn equation. 
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which is the asymmetric dPII equation (|l.2| ). In fact, this hmit was the case that was 
considered originally in the paper [10|. 

However, by itself (3.23) is richer than just d-Pn, a fact that will be made clearer in the 
next section which deals with its Schlesinger transformations. In this subsection, we shall 
examine this equation at the simplest possible level: what are the continuous Painleve 
equations contained in (3.23)? Let us start with the symmetric form of (3.23) which we 
write as: 

Cn+i[x„ + r){rxn + 1) ■ ■ — 7— r - Cn(l - r )x„ 



Xn+l +Xn + r(l + 

XnXn+1 

) Xn-l+Xn + r{l + XnXn-l) 

= Hr{l+2rxn+xl)-r]m{l-xl) (3.38) 

Here and ?7m are the explicit discrete variable related linearly to n respectively m, i.e. 
in general form ^„ = 5n + (^q, r]m = 5m + ryo- The introduction of this scaling of n and 
m is not unavoidable, but it makes the continuous limit procedure more straightforward. 
Next, we introduce the following substitutions: Xn = w, t = ne, 5 = e, (q = rjo = 0, fi = eb, 
Cm = ea, r = e, we obtain at the limit e — > the following equation for w: 

„ w'^ w' 3 2a + 26 + 1 2 2b - 2a - 1 1 
w = Vw w , (3.39) 

w t t t w ^ ' 

i.e. Fill in canonical form. 

We turn now to the asymmetric form of (1.5) which we write as: 

Cn+l{Xn + r){rXn + 1) —77— T " Cn(l " r )x„ r 

Xn+l + X„ + r (1 + XnXn+1 ) Xn-1 + X„ + r (1 + XnXn-l) 

= fir{l + 2rxn + xl) + i/(-l)"+'"(r + 2x„ + rxl) - 7?^(1 - xl) . (3.40) 

Next, we introduce the following substitutions: x„ = — 1 + e'^w for even n, Xn = + 
eH/w + e^u for odd n, where t = e(n — m), 5 = e^, Co = % = li A* = ^^('J — t))/2, 
V = — e^(a + b)/2, r = 1 + e^, we obtain at the limit e — > two equations containing u. 
Eliminating u from these equations we get for w: 

„ w'"^ w' bw^ a 1 /„ 

i.e. Pni although in a slightly non-canonical form. One may remark at this point that 
this continuous limit was obtained for t = e{n — m) and thus in effect, losing one of the 
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parameters of the equation (namely m) . It is not clear whether this is an absolute necessity 
for the continuous limit to exist or not. In the latter case another limit, to some "higher" 
discrete Painleve equation might exist. We shall show in the following sections how one 
can obtain a Py from (1.5) and its Schlesinger transforms. 



3.3 Connection with PVI 



We now show that the similarity reduction also implies a reduction on the level of the 
continuous variables which in our case are the lattice parameters p and q. In fact, from 
the relation (2.28b) in conjunction with the relations given above we can now derive an 
ordinary differential equation for y„ (and consequently Xn) with p as independent variable. 
By construction this equation must be of Painleve type, in fact related to PVI. Starting 
from (2.28b) and its counterpart which can be written as 



^ 1 

p— log Vn,r. 

op 



nan 



mbr, 



(3.42) 



and using Vn+l,m/Vn~l,m — (1 + a„)/(l - a„) together with (|3l2l|) we obtain 

(n + l)a„+i - (n - l)a„__i = m(6„_i - 6„+i) 

'Xn-i+r yn- 



d , (l + an 
-p— log 

dp V 1 - On 



(^-A(-l)' 



i+rn 



nan) 



Vn-i -r Xn+r 

Xn—l Un—l Un 



+m 

Vn-i -r Xn + r 

which by using ( f3.19 ) and the relation = y„,(l — a„)/(l + a„) leads after some 

straightforward algebra to the following first-order equation for a„ in terms of a = a„ and 
y = Vn- 

P{l-r^)^ = (-^ + A(-l)"+™ + na-m){l- a)y + (-/i + A(-l)"+™ + na + m)^ , 
r op y 

(3.43a) 

where we omitted now the label n from all variables. On the other hand, by similar tricks 
we can derive a first-order differential equation for y„, namely by proceeding as follows 
d 

-p— logyn = (n + l)a„+i - nan = + A(-l)"+'" - m6„+i - na„ 



fi + A(-l)"+™ + (^ - A(-l)"+'™ + nan) 



Xn + r 



r yn x^ 
m- 



Xn+r 



-n 



Vn-r 



1 VnXn 

y„(x„ + l/r) y„(l + rx^ 
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where a„ denotes as before the variable a shifted in the m-direction and where we have 
used the counterpart of ( p.21 ) in the m-direction, namely 



rj \yn rj yn 



Thus, we are led eventually to the the following expression for a = an in terms of y = yn- 

{n + m)r(l - y^) - + A(-1)"+™)(1 - r'^)y - p{l - r^)^ 



2na — n + A(— 1 



in+m 



(1 - ry){y - r) 



(3.43b) 



which allows us to derive the following ODE from ( 3.43a ) 



V 



(l-r^) 



d l{n + m)r(l - y^) _ (/, + A(-l)'^+'^)(l - r'^)y - p{l - r2)|| 



dp 



(1 - ry){y - r) 



2n{-fi + A(-l)"+'" + na - m)(l - a)y + 2n(-/i + A(-l)"+'" + na + m 



1 + a 



(3.44) 



where for a we have to insert the expression ( 3.43b| ) in the right-hand side. The nonau- 
tonomous second-order ODE (3.44) for y as a function of the lattice parameter p (which 
notably enters also via r = q/p), is directly related to PVI via 



w{t) = py(t) , i = / 



(3.45) 



in terms of which, taking q = 1 and abbreviating = A(— 1)"^™, we obtain precisely (1.6) 
with the identifications 



a 



1 



1 



{ij, — u + m — n)"^ , (3 = -{jj — u — m + n) 



{^ + 1^ — m — n — 1) 



1 



{li + v-\-m + n+ 1)^ . 



(3.46) 



The connection between the y = yn and PVI (1.6) via ( [j.45D allows us to express 
the general solution of the GDP equation (3.23), or equivalently (1.5), in terms of PVI 
transcendents, namely by writing 



y-a = -Pviip ;a,l3,-i,5;wQ,wi) , 
p 



(3.47) 



and using the identifications (3.46) and the relations ( p. 191 ) and (3.24) to obtain Xn or 
Xn- The initial values wq = w{0) and wi = w'{0) can be identified with initial values for 
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the discrete equation by using (p.43b|) together with (|3^ and §J^). In fact, the latter 



corresponds to a Schlesinger transformation of the form w' = f{w, w, t) , {w' denoting the 
derivative with respect to t and w denoting the shift of the parameters a,f3,'j,6 according 
to n 1-^ n — 1) which ahows us to express for fixed t the initial values of at n = no, 
(given by parameters ao, /3o; TOi (^0 via (3.46)), and at n = no + 1 (given by some different 
parameters ai, 71, (^i) in terms of w and w' at a fixed value of the parameters. Thus, if 
we assume the PVI transcendents to be known functions at each value of the parameters 
for given initial data wq and wi , we have a two-parameter family solutions of the GDP. Of 
course, this starts from the assumption that all two-parameter solutions of the continuous 
PVI are known and fully under control as a class of nonlinear special functions, which is 
obviously not yet the case in this day and age. The initial value problem for PVI was 
considered in from the isomonodromic point of view. 



Remark 3.3: It is well-known that the PVI equation (1.6) also admits a (Lie-point) 
symmetry group G isomorphic to the symmetric group 5*4 on four elements, |^l|, cf. 
also [38|. This group G is generated by three elements gi, i = 1,2, 3, given by 

gi: wh^l/w , t^l/t , {a,P,-f,6) ^ {f3,a,-f,6) 

g2: w^l-w , t^l-t , (a, /3, 7, 5) 1-^ (a, 7, /3, 5) (3.48) 

gs: wh^^w , t^l/t , {a,l3,-f,5) {a,l3,5,-f). 



This symmetry group G was used in e.g. [39, 33], cf. also [|^, to generate one-parameter 
families of special solutions of PVI. 

As can be seen from the identifications (3.46), on the discrete level these symmetries 
of PVI are very natural: the element gi corresponds to the interchange of the discrete 
variables m and n; the composed map gi o g-^ corresponds to the transformation r„^m <-> 
Cn,m) (/^j i^) (~A*i which is the shift associated with lattice parameter pi, = 0. The 
interpretation of the element g2 is slightly more complicated: it can be inferred from the 
combination 

w — 1 

93° 92° 91° 92° 93° 92 ■■ w^t 7 , t^t , {a,(3,-f,6) ^ {6,-f,(3,a) , 

w — t 
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which is equivalent to the replacement y ^ —x and which corresponds to the replacement 
n — > —n — 1. It is well-known, cf. e.g. |2^], that the group of Schlesinger/Backlund 
transformations act on the solutions of PVI by discrete shifts which are the natural shifts 
in the variables n, m as well as n, u. Thus, it seems that the entire symmetry group of 
PVI can be interpreted naturally in terms of operations on the lattice. We will make this 
point more poignantly in the next section where we will derive Schlesinger- and Backlund 
transformations for the GDP. 



Remark 3.4: The bilinear form of PVI was found only recently, [41|, even though r- 
function "forms" were presented earlier, cf. e.g. |2^, We mention here that from 

the bilinear forms in subsection 3.1, namely eqs. ( 3.29| )- (|3.32 ), together with the relations 
( 2.37| ), from which in particular one may deduce 



_ P_( dan,n 

n \ op 



(Tri 



as well as 



n 



^ ' n,r) 



dTn. 



(3.49a) 



(3.49b) 



op op 

and an equation similar to ( |3.49b ) for cr„,m- With the identifications, (taking as before 
g = l), 



w = p 



^n+l,mTn,m+l 



l-(l-p) 



^n,m+l'Tn+l,r. 
W-t p r„_m<7„+l ,m+l 



1 — t I + p an,m+lTn+l,m 

the latter combination corresponding to the action of the element 

w — t t 



(3.50) 



l-t 



t-l 



{a, [3,^,6) ^ (a, (5, 7,/?), 



of the symmetry group G described above, eq. (3.49) together with the total set of 
corresponding bilinear discrete equations for Tn^m and Cn^m, (|3.29D -(^!32|), provides the 
effective bilinear scheme for PVI. In |41| a bilinear system for PVI was obtained consisting 
of eight equations for eight different r-functions, which can be identified with the ones 
appearing (3.50), including equations of second order in the continuous variable t. 
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4 Integrability Aspects 



There are various aspects to the question regarding the integrabihty of the discrete sys- 
tems presented in the previous section. Here we wih discuss basicahy two of them: i) the 
issue of the consistency of the equations from the point of view of initial value problems, 
and ii) the isomonodromy of the lattice Painleve systems. Both aspects pertain to what 
is possibly the most intriguing question, namely: to what extent do these discrete systems 
obey the Painleve property?. The Painleve property being originally formulated for dif- 



ferential equations, [31, 32 1, it has been an important issue in recent years to formulate 



the analogous property for discrete systems, cf. []Tll , 45]. In view of the fact that one 
of the main goals in the original studies by Painleve was to use differential equations as 
a means of introducing new (transcendental) functions with appropriate uniformisability 



properties, cf. [46, 47 [, we believe that the issue of formulating the Painleve property in 
its most general context (namely for difference equations as well as differential equations) 
constitutes one of the key steps in the general problem of developing the proper tools for 
the integration of difference equations. 

4.1 Singularity Confinement and Consistency of the Evolution 

One of the main breakthroughs in the search for the analogue of a Painleve property for 



discrete systems has been the discovery in [11| of the property of singularity confinement 



(SC). In coarse terms SC amounts to the phenomenon of non-propagation of singularities in 
integrable discrete systems. Whilst it is not at all clear whether SC is a sufficient condition 
for integrability, it can be argued that it is at least a necessary condition because it amounts 
to the well-posedness of the discrete evolution of the system. As an algorithmic tool it has 
proven to be very fruitful in predicting the integrability of many discrete systems, notably 
of most of the newly discovered discrete Painleve equations, [12, El]. 



To exemplify the singularity confinement mechanism we apply it to the situation given 
by the constrained lattice system derived in the previous sections. Here, we shall investi- 
gate the nonlinear system by itself without referring to the structures developed in section 
2 from which it was derived. We need as a starting point the general form of the discrete 
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system, allowing for arbitrary coefficients, conditions on which can then be found using 
the SC requirement. In the case of the nonlinear constraint of the form ( f3.17 ), we may 



thus start from the "bare" form of the constraint, |10], 



n ■ ■ \-m — ■ ■ = 0, (4.1) 

'^n+l,m ~l~ '^n—l,m ^n,m,+l ~l~ ^n,m— 1 

i.e. the form that leads to a discrete Painleve equation without free parameters. This 
same form was analyzed in [^] in the framework of the bilinear formalism. However, as 
we now know this is not the most general integrable form of the similarity constraint. Let 



us, therefore, investigate by what we can replace the r.h.s of (4T) to retain integrability, 



taking the point of view that the similarity constraint must be an integrable equation in 



its own right. Thus, taking instead of (4.1) the form 



n ■ ■ h m — ■ ■ = A (4.2) 

where A is in principle a function of n and m. In order to apply the singularity confinement 
criterion we consider the possible singularities of ( [4.2p : either v becomes infinite or it 
becomes 0. In both cases we require that the confinement is in one step, i.e. at the 
iteration immediately following the singular value we already obtain a finite value for v. 
This requirement leads to the following conditions for A: 

[m? — n^) [m{A{n + 1, m) — A{n — 1, m)) + n{A(n, m + 1) — A{n, m — 1))] 
+n[A{n + 1, m)A{n — 1, m) + l){A{n, m + \) — A{n, m — 1)] 
— m[A{n, m + \)A{n, m — 1) + l)[A[n + 1, m) — A{n — 1, m)] 
+2nm[A{n + 1, m) + A[n — 1, m) — A{n, m + 1) — A{n, m — 1)] = , 

and 

(m? — n^) [{A{n + 1, m) — A{n — 1, 'm)){A{n, m + 1) — A{n, m — 1)] 
-\-n[A{n + 1, m) + A{n — 1, m))(A(n, m + 1) — A[n, m — 1)] 
—m[A{n + 1, m) — A{n — 1, m))(A{n, m + 1) + A{n, m — 1)] 
+2nm[A{n + 1, m)A{n — l,m) — A{n, m + l)A{n, m — 1)] = . 

The complete solution of this system is not of particular interest since we have already 
restricted the choice of the coefficients in the l.h.s. of ( |4.2D to exactly m and n while they 
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could a priori also have been any functions of n and m. However, there exists a solution 
that is really interesting in the case at hand, namely a constant A. Moreover, a look 



at (4^) suffices to convince oneself that this constant may have an even-odd dependence 
(since (|4.2| ) involves points that are all shifted by an even number of steps of each other) 
We arrive thus to the precise form given in ( |3.17 ). 

Let us now turn to the issue of the consistency of lattice equation and similarity con- 
straint. For the KdV family where the lattice equations involve vertices around elementary 
plaquettes, whereas the similarity constraints are given in terms of cross-like configura- 
tions, we have thus a system of equations that can be symbolically represented as 



+ 



(4.3) 



As a starting point for the solution of initial value problems for systems of the form (4.3) 
one may take initial data on the lattice according to the following configuration 



(4.4) 



and use the lattice equation and constraint in conjunction to each other to fill the entire 
lattice. The issue of compatibility arises from the fact that at some points on the lattice 
the determination is not unique, and the uniqueness (i.e. single-valuedness of the solution 
around localised configurations) needs to be verified. In fact, this is a discrete manifes- 
tation of the Painleve property, not quite formulated in the same way as the singularity 



confinement phenomenon that was put forward in |11|. Our conjecture is that singularity 
confinement is a consequence of the single-valuedness in the above sense. 

The iteration of the system works as follows: starting with data given on a configuration 
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of lattice points of the form (4.4), one can proceed through the lattice by usinga the lattice 
equation together with the similarity constraint and calculate the values of the variable 
at each lattice point. In Figure 2 are indicated points that are calculated by means of 
the lattice equation by a cross ( x ) and points that are calculated using the similarity 
constraint by an open circle ( o ). However, there are going to be points of possible conflict 
where we can calculate the value in two different ways. The first point where this happens 
is indicated in the figure by a cross within a circle, ( (8) ). For the similarity reduction to 
be consistent, i.e. for the similarity constraint and the lattice equation to be compatible, 
the value at this point needs to be unambiguous i.e. both ways of calculating the value 
must give the same result. 

o X 



Figure 2: Consistency of the constrained lattice system. 



The linear system of partial difference equations 

n (/n+l,m — fn-l,m) + m (/n,m+l — /n,m-l) = A' — A(— 1)""*""* 
(P + Q) (/n,m+l — /n+l,m) = {P — Q) (/n+l,m+l — fn,m) , 



(4.5a) 
(4.5b) 



is possibly the simplest system of the form (4.3) for which one can check the compati- 
bility on the lattice by direct calculation. This system, in fact, is the linearised form of 
the systems of lattice KdV type as given in the previous section. Obvously, it is much 
harder to ascertain the compatibility for the full nonlinear system, and one needs to resort 
to algebraic manipulation programmes (we have used MAPLE and REDUCE ourselves) 
to establish the consistency in a direct way for arbitrary initial data. Thus, by direct 
calculation on the level of the equations, we may state the following: 
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Proposition 4.1: The nonlinear systems of partial difference equations on the two- 



dimensional lattice associated with the KdV family, i.e. either (3. IS) together with (3.14), 
or ^. 1 dj ) together with ( 3.17 ), or ( 3.21 ) together with ([3.2% ) are all compatible, i.e. arbi- 
trary initial data given on configurations of the form (4-4) ^^ac? to a single-valued determi- 
nation of values for the dependent variable on each vertex of the two-dimensional lattice. 



In p^ ] it was argued that the compatibihty of the system should be interpreted as 
a symmetry of the lattice equation. Thus, the lattice Painleve reduction is a source of 
inspiration to a natural formulation of what could be a general definition of a symmetry 
of a partial difference equation. To elaborate a general theory of symmetries of discrete 
equations along these lines is a challenging problem for the future. 



4.2 Isomonodromic Deformation Problems 

In section 2 we have shown how to derive the lattice systems as well as their similarity con- 
straints leading from the infinite matrix formalism. Subsequently, exploiting these lattice 
equations together with the constraints we have shown how to obtain the general discrete 
Painleve equations, either in the form of (1.5) or as (3.23). In this section we will show 
that we can obtain from the infinite matrix structure also the associated linear problems 
which are the relevant isomonodromic deformation problems. In the spirit of the original 



works by Fuchs and Schlesinger [50, 51] the integrability of the discrete Painleve equa- 
tions can thus be assessed in an independent way, namely by exploiting the monodromy 
problems and use them as a starting point for the asymptotic analysis to characterise the 
transcendental solutions, which for the continuous Painleve equations has been an ongoing 
programme during the last two decades, cf. e.g. |l52[| -[^]. A similar programme for the 
discrete Painleve equations has, in our opinion, to be developed as well, but this is beyond 
the scope of the present paper. We contend ourselves here by just providing the tools 
for the further analytic study of these nonlinear difference equations, which at present is 
still only in its initial stages. (Notably, to date the only discrete Painleve equation that 
has been systematically investigated from an analytic point of view is dPI ( |1.1| ), cf. e.g. 
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Mm)- 

We will now show how isomonodromic deformation problems can be derived system- 
atically on the basis of the infinite-matrix structure that we developed in section 2. Ex- 
panding that scheme further we can introduce a set of linear relations for infinite vectors 



which were already defined in ( 2.19 ), namely 



Ufc = pfc (1 - U • ri) • Cfc , (4.6) 



having components n^, {j G Z), with pk given in ( |2.6[ ) and denoting, as before, the 
vector with entries {ck)j = ■ We will use the relations of section 2 to obtain linear 
systems of equations for in terms of the infinite matrices defined there, and then extract 
from these infinite systems finite systems in terms of special entries, leading to concrete 
Lax pairs as well as monodromy problems. Recalling that the action of the operators A 
and I on is given by 

A • Cfc = fccfc , I • Cfc = k-^Ck , 



and making use of the relations ( |2.2| ) we can derive from the matrix Riccati equation ( p.ll|) 
for U the following linear relations for 

{p - k)uk = + A - U • O) • Ufc , (4.7a) 
{p + k)uk = {p-A + V ■0)-Uk . (4.7b) 

From eqs. ( [4.7| ) one can derive Lax pairs for the various lattice equations within the 
KdV family and the relation ( p. 20 ) can, in fact, be used to establish the gauge transfor- 



mations between the various Lax pairs as will be shown below. Next we need to construct 
the similarity constraint for which can be obtained in a similar way, now using the 
similarity constraint (2.17) as a starting point. The result is the following 



^•i. + A 



Ufc . (4.8) 



The next step is to derive concrete Lax pairs and monodromy problems from the 
infinite set of equations encoded in ([4.7D -(4.8). This can be achieved by singling out 
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specific components, the choice of which wih determine which equation within the lattice 
KdV family we are dealing with. In addition we need the relations pertaining to the 
dependence on the lattice parameters p and q. These are given by the basic relations 
consisting of 



d 



1 



1 



op \p + A p — k 



+ U- 



1 



o 



1 



p — *A p + A 



(4.9) 



For the system consisting of the lattice MKdV (3.16) together with the corresponding 
similarity constraint ( |3.17| ) we have the following monodromy problem Q. It consists of 
the linear differential equation in terms of the spectral parameter k, cf. [1lO| ], 

/ 

+ 



-(! + /.) 

7„+, 



2nvn-. 




n—l,m 



+ - 



2mVr, 



1 ~l~ 



I 



n 



n,m+l 







Ipn, 



m—1 ) 



(4.10) 



where 7„ = n + A(— 1)" , in the case of the lattice MKdV, together with the discrete Lax 
pair in the form of the form 



obeying the discrete zero-curvature condition 



(4.11a) 
(4.11b) 

(4.12) 



The matrices L and M are given by 



'l-'n,m{k ) 



fc2 



+ 1 , rn 



(4.13) 



Using eqs. (4.11) to eliminate the backward shifts on V'n.m we obtain an isomonodromic 
deformation problem in the usual sense. In addition, we have the linear differential equa- 
tions with respect of the lattice parameters, namely 



dllJn,r 

dp 



n 



1 



7/. J ihlH 

Yn,m I 

P \ j Vn+l,m + Vn-1,' 



( 



-pVn+l,m , ^ 



n—l,m ) 



(4.14) 
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and a similar equation w.r.t. the parameter q. The infinite matrix scheme of section 2 
provides us with a derivation of the monodromy problems as well. Details of the derivation 
can be found in Appendix A, cf. also The monodromy problem (4.10) was originally 
given for the special values of the parameters = 0, A = 0. The present derivation is also 
taking into account boundary terms which are responsible for the occurrence of the free 



parameters. Our results correspond to the ones of ref. |24|, where it was noted that the 
monodromy problem has the same form as the one for PVI. 

The monodromy problem for the system in terms of Zn^m is obtained from a gauge 
transformation on (4.10), whereas the monodromy problem for the system in terms of 
n„^m, i.e. eq. (|3.27| ) together with ( 3.25 ) is given by the Lax pair 



M 



(4.15a) 
(4.15b) 



with Lax matrices 
/ 



Q — "Un.m+l 



+ * q + Un-. 



(4.16) 



in which the * in each matrix stands for the product of the diagonal entries, in conjunction 
with the differential equation 

/ 



k (pn,m 

dk 



^n+m 

\^ (2A(-1)"+'" - l){np + mq) 7,+ 

/ 



m+l 



+ 



+ 



2np^ 



2p + Un-l,m - Un+1,' 



-1 

yp- Un+l,m. 



2mq 



I 



2q + 



-1 







+ 



''n— l,m 



■'n,m—l 



(4.17) 



\1~ ""n.m+l 

The compatibility of ( [4.11 ) for (j)n,m with the matrices ( |4.16 ) and (4.17) leads to the lattice 
KdV ( 3.25| ) as well as to the similarity constraint ( p.27| ). In addition we have the following 
continuous deformation in terms of the dependence of 4'n,m on the lattice parameters p 
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and (7, namely 





1 



+ 



2np 



op \ 1 n / ' 2p + Un-l,m - "U-n+l.m 




n— l,m 



(4.18) 



and a similar equation in terms of q. 

Finally, a monodromy problem for the SKdV equation, in terms of the variable Zn^m, 
is easily obtained from (4.10) by a simple gauge transformation with the matrix 

1 

l/Vn,n 

using the relations ( |3.15| ). We omit the details here (cf. |^ for the explicit form). 
4.3 Reduced Lax Pairs 



By eliminating the back-shifted vectors Tpn-i^m and ipn,m-i in (4.10) using the Lax repre- 



sentation ( 4.11 ) we arrive at the following linear differential equation with regular singu- 
larities at k"^ = 0,p^,q^ and oo: 



+ 



^n,m ^ Qn.r 



dk"^ \ k'^ k"^ — k"^ — (p- 
in which the matrices 5n,m, -fn,m and Qn,m are given by 



(4.19) 



M+1 
2 



n '"n+l,m^n — l,m _, m ^n,m+l ^n,m— 1 
P ■>}„ + l^m+Vn-l,m Q f n,m+l +f n,m- 1 
1 nv„-l^rn mVn^m-1 



+ 1 ,m ~t~?-'n — 1 ,m ^ri,m+ 1 ^^n,m— 1 



n 



n,m 



Qn,r, 



P{Vn+l,m. 



m 



PVn—l,m 



Q'Vn,m+l Vfi^rn+l'Vn,m—l 



As was noted already in |24] this linear differential equation has the same form as the one 



in the monodromy problem for continuous PVI, cf. |36, 52, 26]. In fact, the continuous 
isomonodromic deformation is provided by the linear differential equations in terms of the 
lattice parameters, namely 



dp"^ 



F 



n,m 



P 

W- 



■ p^ 



(4.20) 
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where 



n 



( 1 



P(^n+l,m ~l~ ^n— l,m) 



V 



and a similar equation for q. From ( [4. 201 ) together with (4.19) we can derive the Lax pair 



for PVI. The discrete isomonodromic deformation is given by the discrete Lax pair ( 4.11 ) 
with ( [4.131 ). Together with (4.19) we have a Lax pair for the constrained system consisting 
of (|3.16| ) and ( |3.17 ). To derive from it an exphcit isomonodromic deformation problem for 
the GDP we perform the gauge transformation 



V'n = V'n("i) 



(4.21) 



1 

which allow us to write the monodromy entirely in terms of the reduced variables Xn, or 
equivalently of ( |3.18| ) which leads to 



K K — t 



+ 



Cn 



Ok \ Hi K — t K — s 

with K = k"^, t = and s = and where the matrices A^, Bn and C„ are given by 

-(^ + 1) n(i_a,) + ^(i_6„) 



(4.22a) 



Ar 







Br, 



n 
2 



( 



l + fln -i(l-a„) 



7„+m - n(\ - an) - m{l - 6„) 
/ 



-p{l + an) 1 



Cn 



m 
~2 



l + 6„ 

-9yn(l + hn) 1 



— (l-6n) 



(4.22b) 



where a„ can be expressed in terms of x„ or ?/„ via eqs. (3.22) and 6„ via the similarity 



constraint mhn = A* — A(— 1)"+'' 



„ . On the one hand we can supplement ( 4.22| ) with 



the discrete equation 



na. 



p 1 



(4.23) 



leading to the eq. (1.5) or with the continuous equation 



dt 



1 



-nan 
n^/t (1 + 







Bn 



K — t 



A. 



(4.24) 



-(n + l)a„+i 

By direct computation the compatibility of the above linear systems can be verified and 
the conlusions can be summarised by the following statement: 
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Proposition 4.2: The system consisting of the differential equation ( 4.2^ together with 
the continuous deformation ( \4-24 ) is compatible iff the ODE (3.44) '^-^ satisfied leading to 



the PVI equation. The same differential equation ( 4-2^ ) together is compatible with the 
discrete deformation ( 4-23i) iff the ordinary difference equation (3.23) is satisfied. Both 
the discrete and the continuous deformation are isomonodromic. 



Remark 4.3: To make the connection with the standard monodromy problem for PVI, 
in the form of refs. |3^, we need to take a different gauge from the one in (4.21), 
namely 

Xn,m = ^n,m'4'n,m j ^n,m — 

^ {Un,m - np - mq)/Vn,m 1 

where the entry in left-lower corner of the matrix Vn,m can be expressed in terms of Vn,m 
using the similarity constraint for Un^m in 

the form (|]2|). Using (Isl^ ) and (2.28a) the 



(4.25) 



continuous isomonodromic deformation ( 4.20 ) then adopts the standard form 



9X71,1 
dt 



K — t 



Xn,: 



Wn,r 



V P v^^ 

'11,1X1^ n,m ^ n m 



(4.26) 



5 Schlesinger Transformations and Self-Duality 



In this section we will derive the Schlesinger and auto-Backlund transformations for the 
GDP equation. This will be linked to the general bilinear scheme for the system. In fact, 
the derivation of the bilinear scheme in this section will be done from a slightly different 
point of view from the one in section 3, leading to a scheme wherein all ingredients of the 
theory will be on equal footing. We must warn that in order to arrive at a self-contained 
description, in this section we deviate slightly from the notations of earlier sections. Thus, 
the variable Xn will not be exactly the same as the one used in ( 3.24| ), but differs by a 
simple fractional linear transformation. 



5.1 A first auto-Backlund relation 



In this subsection, we shall display a first auto-Backlund transformation for the GDP 
which can be obtained in a straightforward way from the basic equations. We start by 
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rewriting them here as: 



— r( ■ ■ j (5.1a) 



Vn+l,m Vn—l.m , Vn^m+l '^n,m—l s 

n h m = fi — u (o.lbj 

"^n+ljin ~l~ 'Wn— l,m '^n,m+l ~l~ '^n,m—l 

where we have assumed that the point (m, n) is such that m + n is even. Next we introduce 
Xn,m = Vn,m/vn+i,m+i and from (5.1a) we obtain Vn,m+i/vn+i,m = {rx + l)/{x + r). Using 
this ratio, we can express (5.1b) (upshifted once in both m and n) in terms of the x's in 
the following way: 

4-1 ■^n+l,m T ■^n,m+l 

+ 1)^SlTI + + 1) r.„.^+i , ^ = M - • (5.2) 

Three more equations can be written centered around the SaniG point ^rj,m' 

3^Ti— l,m~l~^ '^^n,m+l~l~l 



^n— l,m~l~?^ _|_ 'T' '^^n,m+l~l~-^ 

rr 1 -1-1 ~^ •^n^m ™ , , -Ur 



+ ; — +7 \-m— -n = fi + u , (5.3b) 

+ 1 



»^ 3- \- m . = . (5.3c) 

Note that there is some asymmetry in the m, n factors of these equations. This is due to 
the fact that Xn,m is in fact defined at point (n+1/2, m+1/2) whenever (n, m) corresponds 
to the lattice site at which we evaluate fn,m- What is the meaning of Xn,rrJ Clearly, this 
is the same x as the one defined in Section 3, related to the evolution in the n direction. 
The index m here indicates simply the parameter value of the discrete equation and thus 
equations (^^)-( |5l^ ) are just the auto-Backlunds ofthe GDP since they link Xn,m and 
Xn+i,m (or Xn~i,m) to Xn-m+i- These auto-Backlund transformations do not cover all 
the possibilities. In particular, the changes of the parameters v and /i can also intoduce 
Schlesinger transformations. The latter will be examined in the next subsection in the 
framework of the self-dual approach. 



A last interesting feature of the two-dimensional equation ( |5.2| )-( |0| ) is that they can 
be used in order to define one dimensional evolutions in direction different from those of 
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the orthogonal lattice {m,n). Let us consider for example the system of equations (5.3a) 
together with the n-downshifted of equation (5.3b). 



Xn,m~\~'f ^^n — l,m — 1 

ITT ^n— l,m ~Z ; 7Xr~ ■^n— l.m 

-rz ^'^T^ — \ m =fJ'-^, (5.4) 

' '^n^m^ ^ ' J-n— l,m — IT^' ' 

(notice the change of sign in the r.h.s.). 

Introducing the one-dimensional notation X„ = Xn,n+k, Yn = Xn-i,n+k (with k = 
niQ — no, where the index characterizes the starting point), we obtain thus what we call 
an asymmetric form (the name has its origin at the asymmetric family of QRT mappings 



introduced in ||6^, cf. also pV\ 



(Cn-Ac) ^+7— — + (Cn+i + K)7^^7+^-— =/i + i^ , (5.5a) 



X„+r rX„_i+l y 

rX„+l , , X„-i+r -^'^ 

with again ^„ = (5n + Co 2k = 5k. This mapping is a discrete form of Py equation as 
can be assessed through the continuous limits. Indeed, putting X = 1 + e{w + 1)/(2 — 2w), 
Y = 1 + eu, 6 = e, r = —{1 + e/A), fi = eg, Co = 0, u = eh we obtain: 

„ ,2/ 1 , -L N ^' , jw - 1)% , b cw dw{w + 1) 

«^ ^ + 7 + 2 01" + -) + — + — , (5.6) 

2w w — 1 z w z w — 1 

i.e. Pv with z = ne,a = {k + g-hf/S, h = -{k- g + hf/8, c = {g + h)/2 and d = -1/8. 
At this point it is interesting to show that the equations by ( ^.5[ ) are equivalent to 



another discrete form of Py, obtained in |61|. We start by introducing the new variables 
through the following change: 

y„ = + ^ , (5.7a) 

V T 2En + 6 + 2K- fJ,-U 

Xn = -l^n+l-Z= 77 ^7 — • (5.7b) 

2r.n - AC - 36 - 2k - - u 



Prom equations ( p-5\ j we obtain, in terms of the variables H„, T„: 

_ (2H„. -4Cn-36 -2k- u){2En - 4Cn - S + 2k + fi + u) 

" (2H„ + 5 + 2K-fi- u){2En -6-2K + fi + i^) ^^'^''^ 

„ _ 2Cn + ^l-l^ 2Cn + - . 

"--^ + - T„/r + 1 + T. + r " ^^"^^^ 
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Introducing the quantities a = fi — v, h = k + [5 — ^ — v) /2 and c = k + {5 + ^ + u) /2 we 
can write the equations as: 

i n i n+1 = 7= , (b-Uaj 



This system coincides precisely with the one obtained in |£T| from the degeneration of 
g-Pvi- It goes without saying that the continuous hmit of (^.90 is again Py- 



5.2 Bilinear expressions and Hirota-Miwa equations 



The system of equations ( |5.2| )-( |OD offers a fully symmetric description of the evolution 
along the direction n and m. However, the GDP has clearly more Schlesinger transfor- 
mations (those corresponding to the parameters /i, v) which do not seem to fit into this 
symmetric scheme. We shall show in this subsection that it is possible to introduce a fully 
self-dual description of the GDP in which all four directions n, m, ^, v play the same role. 
The key to this approach is the bilinear formalism. Let us start by the bilinearization 
of equations ( |5.2p and ( |5.3| ). We will find it useful to distinguish between four different 
r-functions will be necessary: F, G that exist only for 'even' lattice points (i.e. n + m 
even) and H, K present at the 'odd' points, namely 

Vn,m = and Vn^\,m = • (5.10) 



Thus, comparing with (|230| ) the F and H can be identified with the r+ at those sites and 



the G and K with the t_ (or, equivalently, the an,m- Next, we substitute into ( 3.16 ) and 
obtain: 

Hn,m+l-l^n+l,m{pPn,mGn-i-l,m+l ~l" QPn+l,m+lGn,in) 

= Hn+l,mKn,m+liQFn,mGn+l,m+l + pFn+l,rri+lGn,m) (^-H) 

For this expression we see that we can separate provided we add the same quantity 
— (p -|- q)[Hn+i,mKn,m+iHn,m+iKn+i,m) to both sides of the equation. (The choice of 
the multiplicative factor is based on the convention of Hirota who takes the sum of the 
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coefficients of the Hirota-Miwa equation to be equal to zero). Separating (5.11) we find. 

pFn,mGn+l,m+l + qFn+l,m+lGn,m = {P + q) Hn+l,mKn,m+l (5.12a) 
pFn,TnGn-l,m-l + QFn-l,m-lGn,m = (P + Q)Hn-l,mKn,m-l (5.12b) 

where (5.12b) is in fact the down-shift in both m and n of what we obtain in separating 
(5.11). In eq. (5.12a) we recognise immediately (3.31a). 

Thus, the autonomous mKdV equation is bilinearized to a system of Hirota-Miwa 
equations. A close look at (5.12a) and (5.12b) reveals that the two equations, although 
non strictly identical, describe the same evolution. Simply the direction of evolution is 
reversed between (5.12a) and (5.12b). 

Next, we turn to the bilinearization of ( 3.17] ). Since the equation relates f's at 'odd' 



points of the lattice (shifted by one with respect to Vn,m) only the H, K functions enter. 
We find: 

Hn+l,mKn~l,m{{n + m — fl + u)Hn,m+lKn,m~l + {n — 171 — fl + u)Hn,m-lKn,m+l) 
= Hn-l,m.Kn+l,m{{n - 171 + fl - u)Hn,m+lKn,m--l + {n + 171 + fl - v)Hn,m-lKn,m+l) 

(5.13) 

Again we separate by adding the quantity gHn~i,mKn+i,mHn+i,mKn-i,m to both sides of 
the equation. This time it is not convenient to implement the Hirota convention and thus 
we just choose g = —2n. 

{n + m- fl + u)Hn,m+lKn,m-l + {n - 171 - fl + u)Hn,m-lKn,m+l - '^'riHn-l,mKn+l,m = 

(5.14a) 

{n-m + fl - u)Hn,m+lKn,m-l + {n + m + fl - u)Hn,m-lKn,m+l " '2nHn+l,mKn-l,m = 

(5.14b) 

(Notice that both equation are invariant with respect to the parity n —n). Again 
the two equations, although not identical, describe the same evolution with the direction 
reversed. Expressions (5.14a) and (5.14b) are reductions of a more general Miwa equa- 
tion which involves all four products Hn+i,mKn~i,m, Hn-i,mKn+i,m, Hn,m-iKn,m+i and 
Hn,m+iKn,m-i and has a free coefficient (which allows the reduction to be constructed). 
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Since we already presented the bilinear form of GDP in section 3 it makes sense to 
compare it to the one obtained here. We have used various notations for the r-functions 
so let us first establish the equivalence between them. Comparing (|2.30| ) as it was reinter- 
preted in Section 3.1 with ( 5.10 ) we have: Tn^m = Gn,in-, (^n,m = -^n,m for n + m even, and 



Hn.m for n + m odd. 



It is straightforward to show that ( ^.14 ) comes directly from ( 3.32| ) provided we express 
o'n,mTn,m using ( 3.30 ). Similarly ( ^.12 ) can be obtained from (|3.29|) -( 3^ ) after some 
algebraic manipulations. Conversely, we can obtain (|3.29D - (|^) starting from (|1|) 
and ( ^.14 ). We start by noting that both sets of equations are invariant under the gauge 
r — > ec(n2+m2)^ where r stands for any of the F,G,H,K. On the other hand the l.h.s. 
and the r.h.s. of (|3.30| ) behave differently under this gauge. Thus, we can use it in order 
to ensure that (3.30a), for instance, is satisfied at a given point {n,m). Once (3.30a) is 
satisfied at one point one can use ( ^l2D and (|1|) to prove that (|]2|), (|3^ ) and ( ^) 
hold true at each point of the two-dimensional lattice. 



5.3 Self-dual formulation of the GDP 



In the previous subsection we have seen how the GDP can be expressed as a system 
of bilinear autonomous and non- autonomous Hirota-Miwa equations. Although offering a 
complete description of the system at hand these equations treat in a completely asymmet- 
rical way the dynamic variables (n, m) on the one hand and (/i, u) on the other. However 
the latter play precisely the same role as the former since their evolution is also related 
to Schlesinger transformation. What is needed is a symmetrical, self-dual treatment of all 



four parameters, in the spirit of the "Grand Scheme" of |27|. 

As a first step towards this self-dual description, we introduce a change in the notations. 
The variables n, m, /u and will in what follows be represented by ni, n2, n.3 and 714. Next, 
instead of the four r-functions F, G, H and K, we introduce a single function r which 
is a function of all four variables, r = T(ni, n2, ?t-3, "^4). Finally we use the short-hand 
notation introduced in [ p7| , where upper and lower indices are used in order to denote 
shifts r'^ = r(. . . , rifc + 1, . . .), = t(. . . , — 1, . . .), (with the obvious disadvantage that 
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a shift of j steps requires j repetitions of the index of the respective variables). 

Prom the results of the previous subsection 5.2, cf. also subsection 2.3, we know that 
f is a ratio of two r's and it turns out that the good choice is to consider only r-functions 
with an even number of shifts. We have thus the following identification of the r functions 
of the previous subsection (it goes without saying that there is some arbitrariness in the 
choice of this correspondence): 

Fn.m = T ) Gn,m = "T-, -f^n+l,m = ) I^n+l,m = T ■ (5.15) 

Note that for the r-functions F and G the sum of the shifts in the directions 1 or 2 is 
even, while this sum is odd for H and K. We can now transcribe equations ( ^.12 ) to: 

p^34^12 ^ ^^1234^ ^ ^ ^)^13^24 ^ (^^^g^) 

pT^\i2 + qrllr = {p + g)rf t| (5. 16b) 

Similarly, the non-autonomous equations ( |5.14| ) are written as: 

(ni + n2 - ns + n4)T^^r| + (ni - n2 - ns + ni)T2T'^^ - 2niTfT^^ = (5.17a) 
(ni - n2 + ris - n4)T'^^T2 + (ni + n2 + ns - n4)T^T^^ - IriiT^^rf = (5.17b) 

This system is not sufficient in order to describe the evolution in all four directions. We 
can in principle extend these particular forms formally to the remaining dimensions by 
permuting the indices, but this leads to incompatibilities. So we must first modify ( ^.16| ) 
and ( 5.17| ), introducing the appropriate gauge factor r (pr. To make a (very) long story 



short we give here the factor (p 

ai + 04 / V «2 + 04 / V as + oi / V as + 02 / V as + 04 / 



n^+n| n'^+n'^ n'^+n'^ n'^+n^ 

,2 ^2\ — 4 / „2 „2\ — 4 /„2 ^2' ' " 



ag - gj \ / as - a2 \ / a^ - 04 \ / aa - 04 

aios / V a2a3 / I 0104 / I 0204 



with: 



ai - 02 V"! ~ "4V°'3 ~ "2 + V"2 ~ '^4V"3 X , , 

SZ — ; ; ; ; ; , (5.19) 

ai + «2 , /„2 - ' " - ' - - ' - 



'«! - «4V "3 - «2 - V "2 - «4V "3 - «1 
where the Oj, i = 1,2, 3, 4, are any four numbers subject only to the constraint 



p \i CLo a, A 



(5.20) 
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The careful application of this transformation leads to the system consisting of: 

(ai - 02) (03 - a4)r^V^'' + (ai - 03) (04 - a2)T^V^'' + (ai - 04) (02 - a3)r^^r^^ = , 

(5.21a) 

2nia2(ai - a3){ai - a4)ri3r^'* + (ni + 722 + 113 + 724)01(02 + a4)(a2 + a^)T2'^l 

-{ni - n2 + n3 + 724)01(02 - a4)(o2 - 03)r23r^'* = , (5.21b) 

which must be understood not only as they stand, but also with all possible permutations 
of {1,2,3,4} and also any number of the parities — > r^, o^ —o,k, — > — ti^. 
Equations (5.21a) and (5.21b) (augmented as we just explained) are obviously self-dual. 
All four directions play the same role. Thus, the GDP can be understood as defining an 
evolution in a four-dimensional discrete space Z^. 

5.4 The remaining auto-Backlund transformations 

In order to introduce the auto-Backlund transformations for the remaining directions, 
(/2, i^) we start with some considerations of singularity structure. Using the bilinear for- 
malism introduced in the preceding subsection, we have: 

^ = :j;t2 = 771234 (^•^^) 



The second special expression involving x is also simply given as 

(5.23) 



X + r t^^t'^^ 



ra; -j- 1 v"^ -^14^23 

where r = q/p. From these two expressions, one concludes immediately that the four 
singular values of x are 0, 00, —r and — 1/r. As we have seen in subsection 4.1 the auto- 
Backlund transformations between directions 1 and 2 involve directly expressions ( p.22| ) 



and (5.23). Let us now consider the auto-Backlund transformations that relate directions 



1 and 3, i.e. the Schlesinger transformation for 723 while the evolution is still along rii. 



Permuting 2 and 3 in expressions ( 5.22 ) we find that, up to a gauge, the relevant quantity 
is: 



^13^24 



y = 771234 (5-24) 
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One can understand the origin of this gauge factor in the fohowing way. When starting 
from the fully self-dual equations, one must introduce a different gauge whether one wants 
to specify an evolution in the couple of directions (1,2) or (1,3). Thus the relative gauge 
of y with respect to x comes from the combination of these two gauges. Comparing ( ^.241 ) 
and (|]2|), and using ( |5.16[ ), we conclude that y oc (x + r) with a multiplicative factor that 
remains to be fixed by the appropriate gauge. Instead of computing this factor by carefully 
going through the gauge transformations we use a short-cut. We start by remarking that 
the permutation of indices corresponds to the permutation of the singularities. Let us 
give an example. As we have seen, the permutation of indices 2 and 3 transforms x to 
x + r (up to a global factor). Thus, the singularities {0, oo, — r, —1/r} are transformed into 
{r, oo, 0, r — 1/r} but we know that, in the appropriate gauge, the product of the two finite 



singularities must be unity. So we must renormalize them by a factor A = l/Vr^ — 1, 
which leads to y = {x + r)/\/r'^ — 1 having its singularities at {— s, oo, 0, — 1/s} with 
s = —r/\/r'^ — 1. 

We can now establish the auto-Backlund transformations involving direction 3. We 
find: 

+ ' y::::s . + (/^ + ds^--^ = - - - > (^-^s) 



which is the equivalent of equation ( |5.2| ) of subsection 4.1. The analogs of equations 
( |5.3| ) can be established in a straightforward way. The auto-Backlund relations involving 
direction 4 can be obtained in a similar way. 



6 Conclusions 

In this paper we have established a direct connection between the fully discretised lattice 
systems associated with the KdV family and Painleve equations via the similarity reduction 
on the lattice. In particular, under the reduction PVI emerges as the continuous equation 
in terms of the lattice parameters. It is well-known that Schlesinger transformations act 
on the solutions of PVI in a discrete way , cf. e.g. [^], and it is therefore not surprising 
that a lattice structure sits at the background of this equation. However, by making the 
connection with the lattice KdV we have been able to identify this lattice structure in 
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explicit terms, and in its most attractive form it leads to an explicit difference equation, 
the GDP (1.5). This equation, as an integrable difference equation of "Painleve type", 
can then be studied in its own right. Furthermore, the relevant bilinear structure of PVI, 
found in [^ ], seems to fit naturally in the discrete structure demonstrating that, in fact, 
the continuous and discrete aspects for these systems are intimately interlaced. The self- 
dual bilinear structure of the discrete equations exhibited in subsection 5.3 demonstrates 
that independent variables and parameters of the Painleve equations live on equal footing. 

We expect that the results highlighted in this paper will allow to gain new insights 
in the nature of both the discrete as well as the continuous transcendents. The rather 
complicated combinatorics that is revealed by some of the derivations, e.g. in subsection 
3.1 and in the patterns of section 4.1, we believe are indicative of the existence of some 
highly nontrivial addition formulae for the transcendental solutions, reminiscent of the 
contiguous relations for functions of hyper geometric type. Ultimately, we hope that these 
results will form an inspiration for a future analytic theory of ordinary difference equations. 
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Appendix A: Derivation of Monodromy Problems 



We present here some of the details in the derivation of the monodromy problems for the 
GDP, (1.5), which was also given in the form (3.23). We will derive the relevant formulae 
from the infinite matrix structure of section 2. All Lax pairs and monodromy problems 
can be derived from that structure. 



a) Lattice (potential) KdV: In this case we take the components and to 
form the two-component vector 



ct>k^{p- knq - ky 



(A.l) 



for which it is immediate to derive the linear relation 

p — u 



L 



k(Pk 



p + * p + u 



(A.2) 



in which u = Uo,o as before and * stands for the product of the diagonal entries of the 
matrix Lfc. Eq. (A.2) is one part of the Lax pair of the potential lattice KdV equation. 
The other part is identical apart from the replacements: u u, p ^ q, leading to 

/ 



4>k = Mk(t)k , Mfe 



1 



+ * q + u 



(A.3) 



In this way we obtain the Lax representation for ( 3.25| ) from the compatibility condition 

In order to obtain the monodromy problem that gives us the similarity constraint we have 
to use (4.8) and write down what this gives us for the vector cpk- Applying the constraint 
for and uj. respectively we obtain 



d 



1 



ul - v-pu-^^ I + mq 



k^,4 = K-ir-uUnp^^_^ 



k-ul = (1 - A(-l)"+-) nl + 2A(-l)"+-(pn + mq)ul 



1 



q — k 



7/0 7, 7/^'*') 



q"-k 



(A.4a) 



+np 



,1 „,o 



p — k 



ip) 



Uk-Uk- s-pul ' ]+mq 



q — k 



(A.4b) 
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where u^^^ is defined by 



u 



a + A / 



(A.5) 



and where the objects Va and Sq, were defined in ( |2.21b| ). In order to bring this set 
of equations into shape we need a number of relations, namely ( p.22| ) and ( 2.25| ) for 
a = —p, —q together with (|2.39|) and ( p.40|) , as well as from (4J) the relations 



(a) 



(A.6) 



applied to a = p, (and similarly to the equations with q). 

For the dependance on the Backlund parameters p and q we have the following relations 



d_ 

dp 
d_ 
dp 



,{p) 



1 



Ui. = n \ Uu + s^„u 



p-k " 
(p) 1 



(A.7a) 
(A.7b) 



p — k 

with similar relations describing the dependence on q. Thus, we obtain a linear system of 
the following type 



(Ai 



dp 









4>k+n 


I 

V-pVp 




n 






:) 










^ s—p'^p 





and similarly a linear system in terms of q. 

The result is the following differential equation for the vector (pk 
( 



n + m 







^ —np — mq n + m + 1 
/ 




2p + u —u 




2p + u — u 



(A.9) 



which is (4.17. We must next use the inverse of the Lax representation (4.11), ( [4.15| ) to 
express (pk and (j)k in terms of cpk , thus yielding a pure linear differential equation in terms 
of the spectral parameter /c^. 
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b) Lattice (potential) MKdV: To obtain the Lax representation for the lattice MKdV 
equation a similar calculation need to be performed. In this case we single out the com- 
ponents together with for some fixed value of a, to construct the two-component 



vector 



i^k = {p- kT{q - ky 



(A.IO) 



V 



First, using the relations ( p. 22 ) we obtain from (4.8) now the Lax matrices in the form 



p — a 



Va 



(A.ll) 



and a similar expression for by making the usual replacements. 



Remark A.l: We may use eqs. ( 2.22 ) to show that the Lax representations ( [A.2[) and 
( A. lip are related via the following gauge transformation, 



L 



{MKDV) 



UkL 



{KDV)jj^l 



See '^a 

k^-a^ 



(A.12) 



(and a similar formula relating the matrices M^). 



Now, the dependence on the Backlund parameters p and q of the quantity u^"^ plays 
a role, and this is given by the relation 



dp 



u 



(a) 



n 



a 



1 



U 



(P) 



(A.13) 



. (p — a)(p — /c) Y X- p — a 

with a similar relation describing the dependence on q. Thus, we get for the vector ipk for 
Q = the following linear system 



dp 





^ 1 ' 


\ 1 


n 










/ V 



-(1+pSo,-; 

pv-pVp 



(A.14) 



and an accompanying linear system in terms of the dependence on q. 

In order to derive a monodromy problem from (4.8) we need to restrict ourselves to 
the case a = 0, because otherwise we would obtain from the term with the operator I 
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derivatives with respect to a, hence equations which are no longer closed-term. Using the 
following additional relation for u'^'^ = u"i^~^^ 

fc|4°^ = -(^ + 1)4°) + n (-A^4°) + (1 +p.o,-p)4^^) 



+m I -^uf + (1 + gso,-j4'^ 



(A.15) 



in addition to (A. 4a) and the relations (2.39a) and (2.39c) we easily obtain the following 
monodromy problem 





/ 



\ (00 

\ i^k + {n + m) \ I V'fc 

A(-l)"+'" ) 10 1 



+■ 



vq + Vq 



Wo 1 , 2mvQ 

i/Jk + 

-p I ~ vo + vo 



( 



V 







(A.16) 



which is (4.10). Again, in order to obtain the purely differential equation we need to use 



the inverse of the discrete Lax pair ( A. 11 ), and its counterpart, to replace the backward- 
shifted ipk and ipk- 
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